From factorizable S-matrices to Conformal 

Invariance ^ 

Roland Koberle 
Institute) de Ffsica e Quimica de Sao Carlos 
Universidade de Sao Paulo 
Sao Carlos 13560 
BRASIL 

February 1, 2008 



Tectures delivered at the Jorge Andre Swieca Summer School, 1990. 
t Supported in part by CNPq. 







1 Introduction 



Recently there has been spectacular progress in 1+1 dimensional quantum 
field theory 0]. Although one may argue about whether this field belongs 
to physics or mathematics, it is abundantly clear that a successful non- 
perturbative attack on 4-dimensional field theories will require a thorough 
understanding of the simpler counterpart in two dimensions. This is about 
all we will say to justify our living in so narrow an environment, except per- 
haps to proclaim our believe that if someone finds beautiful structures in 
mathematical physics, we may be sure that Nature will make use of it. And 
beautiful structures surely we will find. In the present lectures I'll try to 
make them as accessible as I can. 

Trying to make progress in physics requires in general hard work. Yet this 
can be very rewarding, if miracles do happen along the way - and in exactly 
soluble models they do. After all some fine tuning of coupling constants 
has taken place in order to guarantee solubility. One way to ensure this is 
to impose factorizability on the S'-matrix, i.e. we require the n-particle S- 
matrix to be a product of 2-particle ones, consequently not allowing particle 
production. This is guaranteed, if the model possesses an infinite number of 
conservation laws 0, DJ. 

In section we will characterize our models by specifying their conserva- 
tion laws and then explain the Bootstrap Principle to construct factorized S- 
matrices. This program is executed in section |3] for some models. In section ^ 
the Thermodynamic Bethe Ansatz will be used to connect these S"-matrices 
to their ultra-violet limiting conformally invariant field theory. Section [5] con- 
tains a short review of conformal invariance, so that the unfamiliar reader 
doesn't have to go shopping somewhere else, although in practice he may 
want to look up ref . |J . Finally in section ^| we will study perturbations of 
conformally invariant field theories, find out which of them are candidates 
for exhibiting an infinite number of conservation laws and thus make contact 
with section For the benefit of the reader - hopefully -, I have tried to 
make these notes as self-contained as possible, at the risk of copying from 
too many places already well known material. 

Prom a logical point of view, this last section should come first, but it has to contain 
older, more technical material for which excellent review articles exist 0, Q[. Therefore I 
have placed it near the end, hoping to delay as much as possible the loss of my audience. 
In this spirit I tried to exhibit - as we say in Portuguese - 'o caminho das pedras'. 
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2 1 + 1 dimensional factorizable 5-matrix the- 



ory 

2.1 Kinematics and scattering 

We will consider scattering of n particles A a , a = 1,2 ... ,n, whose masses 
are m a . Their momenta satisfy the mass-shell condition 

P^ = PP = m 2 , (1) 

where the light-cone components of p^ are p = p° + p 1 and p = p° — p 1 . 
The on-shell condition may be conveniently parametrized introducing the 
rapidity 9 : 

(jAp 1 ) = (m cosh m sinh (2) 

In order to ensure exact integrability we will assume that the field theory 
giving rise to our ^-matrix possess an infinite number of nontrivial, commu- 
tative integrals of motion. 

Although this is a well suited requirement in d = 2 dimensions, for d > 2 
a theorem by Coleman and Mandula states that in a non-trivial Poincare 
invariant field theory the most general invariance group is a product of the 
Poincare and an internal symmetry group. The proof assumes some analyt- 
icity in energy and momentum transfer ( angle ). This poses no problem in 
two dimensions, since the scattering angle is anyhow only either or n. 

These non-trivial, i.e. other than energy-momentum, conserved charges 
Qm,...,ij, s transform as s-th order tensors under the Lorentz group. In the 
light-cone representation we call them P s ,s = si,S2, . . ., where the label s 
indicates the spin of P s . Actually in two dimension we have no rotation 
group, so that spin refers to Lorentz- spin, specifying how P s behaves under 
a Lorentz transformations L a : 9 — > 9' = 9 + a . Thus 

P s - P' s = e sa P s . (3) 

For example the momentum p has spin one : p — >• p' = e a p and the parity 
transformed p has spin minus one. Therefore we have P\ — p and P_i = p. 
Since parity relates the integrals of motion P„ s to P s , we will consider only 
s > as we will deal only with parity conserving theories. 
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P s acts on one-particle states as 



P s \ A a (p) >= uj a s (p) \ A a (p) > . (4) 

Since P s carries spin s the Lorentz-transformation property equ.( |3[) requires 
u}g(p) to be of the form 

<(p) = <P S = <{rn a ) s e s \ (5) 



or 

P s \A a (p) >= ^\A a (p) >, r s = <«T, (6) 

where K a s are real constants. For example nf = 1. 

P s are integrals of local densities. Therefore their action on well separated 
multiparticle in or out states is the sum of the one-particle contributions. 
Consequently in a scattering process 

p u ...,p n -► pi,...,p^ (7) 

we have 

n m 

Efe) s = E(^) s - (8) 

8=1 1=1 

At least one non-trivial conservation law with s > 1 implies in the above equa- 
tion n = m, ° thus forbidding particle production and only time-delays and 
exchange of quantum numbers are allowed 0. Therefore after an eventual 
relabeling, we have p$ = p\ and a product of energy-momentum conserving 
(^-functions can be factored out of both terms of the r.h.s. of S = I + iT . 
When writing S a b(s) we will always assume these 5-functions to have been 
factored out. 

The n-particle S'-matrix is then a product of n(n — l)/2 two-particle S- 
matrices. This decomposition can be effected in several ways and all of them 
must give the same result. The ensuing consistency conditions are called 
Yang-Baxter factorization equations |]. [J. 

Let us order our particles along the line as p\ > P2 > ■ ■ ■ > p n , so that 
in-states are arranged along the space-axis according to decreasing momenta 
and out-states the opposite. Thus m-states are going to scatter, whereas 

2 Except at most for a discrete set of momenta, which we exclude appealing to analyticity 
of the S'-matrix 0. 
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particles in out-states are running away from each other. The action of the 
S-matrix on a two-particle state is then defined as 



\a{9 1 ),b(9 2 ) >= 

S:' b h '(9 1 - 9 2 ) | a'(9 2 ), b'{9 1 ) > + S b ^'(9 1 - 9 2 ) | b'(9 2 ), a' {6,) > . (9) 

Note that the final states are ordered opposite to the initial states in rapidity 
space, so that the S'-matrix carries m-states into owt-states. This ordering 
introduces at most a phase, which is irrelevant for two-body scattering, but 
when we write a many-body amplitude as a product of two-body amplitudes, 
this phase convention simplifies the factorization equations. 

For simplicity we assume that in eq.(|9]) the second reflection term van- 
ishes and only the first transmission term is present. Unitarity requires 
then the S'-matrix element = S a b to be equal to expzS(9) in the physical 
region. The Yang-Baxter factorization equations are now trivially satisfied 
for these diagonal S-matrices. One way to ensure vanishing reflection is to 
require, that all particles A a have different masses, whenever a ^ b, i.e. there 
is no degeneracy and no internal symmetry present, ( in particular particle 
and anti-particle are identical A a = A a ). In this case no exchange of quantum 
numbers is possible. Or, if this is not the case, some other mechanism, such 
as antiparticles being bound-states of particles [[J, particular conservation 
laws etc. has to be invoked to ensure vanishing reflection. 

In order to pin down the possible functional dependence of S(9), we have 
to discuss its analytical structure, which turns out to be very simple. In terms 
of the Mandelstam variable s = {pi+p 2 ) 2 S(s) has cuts along s > (mi+m 2 ) 2 
and s < (mi— m 2 ) 2 , which are required by two-particle unitarity. We suppose 
these to be the only cuts in the complex s-plane, since there is no production 
and anomalous thresholds turn into poles in 1 + 1 dimensions [[J. The cuts 
are eliminated by the uniformizing variable 9 : 

s = s(9\ — 9 2 ) = ml + m\ + 2m\m 2 cosh #12, (10) 

where 9\ 2 = \9\ — 9 2 \. The mapping from the s- to the #-plane is shown in 
figure pTT| . 

Since poles in the s-plane generate poles in the 6*-plane, S a b(9) is a mero- 
morphic function of 9. 

Real analyticity in the s-plane ( Schwartz reflection principle ), S*(s*) = 
S(s) , becomes S*(— 6*) = S(9) in the ^-plane. In particular we see that 
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the bound-state region (mi — m 2 ) 2 < s < (mi + m 2 ) 2 , where S(s) is real, 
is mapped into the segment < 3m0 < 7r of the imaginary #-axis. The 
unitarity condition <S(s)<S f (s) = &(s)S(s) = I or S(9)S\9) = S^(9)S(9) = 
J, then becomes for real 9 



s ab (e)s ba (-o) = 1. (ii) 



Parity invariance imposes 



S ab (0) = S ba (9) (12) 
and charge-conjugation invariance implies 

S ab (0) = S a - b (6), (13) 

where b indicates the antiparticle of b. Finally crossing-symmetry 

SM = S m (s) = S ab (2mj + 2m 2 - s) (14) 

becomes 

Sa b {6) = S aS (6) = S ab (m-6). (15) 

It follows from equs. (^) and (0), that S(9) is a 2m- periodic function of 
9. We now exhibit the general solution S ab (9) of these equations. Since we 
know that the scattering amplitudes are bounded functions of the momenta, 
one can show that this implies that any meromorphic, real analytic, 2iti- 
periodic function f(9) satisfying f(9)f(—9) = 1, can be represented as 

aeA sinh[l/2(0 - tira)] 

where A is a set of complex numbers invariant under complex conjugation. 
If there are no unstable particles present, then all poles of f{9) occur on 
the imaginary #-axis, implying a to be real and we may restrict it to —1 < 
a. < +1. In this case f(9) has a simple pole of residue 2«sina7r at 9 = iair 
and a simple zero at 9 = —lair. Besides this we note the following useful 
properties: 

/«(*) = fa+2(9) = f- a {~9) 
fa(9)f- a (9) = U9)U-9) = 1 

/a(«T-0) = -/l-a(0) (17) 
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Figure 1: Mapping from s- to #-plane, showing the position of cuts and 
bound-state poles. 
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Figure 2: Diagrams responsible for poles in the direct and crossed channel. 

f a (0 ~ ™P)U0 + m(5) = f a+p {9)f a _ p {9) 

f (9) = -h{9) = 1 . 

If at least one of the particles A a or A b is self-conjugate, then crossing sym- 
metry 

S ab (9) = S ab (m - 9) (18) 
implies, that up to a sign S a b(9) must be of the form 

p (q\ — t (f) \ f m _ sinhO + z sin an _ tanh[l/2(g + tair)} 

sinh9 — ismair tanh[l/2(6' — ia%)\ 

The functions F a (9) satisfy: 

F a {9) = F a+2 (9) = Fi_ a (9) = F_ a {9) 

F a (9)F_ a (9) = 1 (20) 
F a (9 + mf3)F a {9 - mf3) = F a+p {9)F a _ p {9) 
F (9) = 1 . 

For < a < 1/2 , F a (9) has simple poles at lan and i(l — a)n with residues 
2«tan an and — 2«tan an, respectively, as well as zeroes at —tan and — i(l — 
a)n. F\{9) has a double pole at in/2 and a double zero at —in/2. 

Since bound-states will be all important in setting up the bootstrap prin- 
ciple, let us have a fast look at them, both in the s- and the #-plane. In the 
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s-plane S ab (s) has bound-state poles for < s = m c < (m a + mj) , whose 
residues are positive due to unitarity : 

S*(*)~(&)7(*-w#- (2i) 

When we consider S a b{@), we suppose always that the (^-functions of the form 
Y[S(6i — 6j) have been factored out. Therefore in going from S(s) to S(6), 
we have to take into account the Jacobian 

^ n a,P n b } = m a m b smh(8 ab ). (22) 
0{t> a , Vb) 

Therefore near the position of the pole in the #-plane given by 

m 2 c = s (6 = iu c ab ) =m 2 a + m 2 b + 2m a m b cos u c ab , (23) 

i.e. near 6 = iu c ab1 S ab {6) behaves as 

S ab (6)~ t Rl b /(8-tu c ab ), (24) 

where R ba = f^bfab- F rom crossing-symmetry we see, that there is also a 
pole with negative residue — R c ab at 6 = m c ab = i(n — u c ab ). This behavior is 
caused by the diagrams displayed in figure 0, which shows, that f abc = f£ b 
is a completely symmetric function of its indices. From the law of cosines 
applied to the triangle with sides m a , m b , m c , we obtain from equ.(p3|) 

<b + 4a« + <o = 2tt. (25) 

Up to now we have listed the properties of S ab (8) following from general 
principles. These are of course insufficient to construct the two-particle S- 
matrices. Since the factorization equations are void for diagonal S'-matrices, 
we need another input. This is furnished by the Bootstrap Principle . Sup- 
pose our theory contains the set of particles Ai, . . . ,A n . Then we require 
that the bound-states of S ab (9) with positive residue represent one of the 
particles A a and vice-versa - any particle A a occurs as a bound-state in some 
S- matrix element. 

Suppose a particle A c occurs as a bound-state in the scattering of A a A b 
at 6 ab = iu c ab . Then scattering A c with A d must give the same result as the 
three-particle scattering 

Sabd(0 a ,8 b ,6 c ) = S ab (9 ab )S a d(9ad)S b d(0 b d) (26) 
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at the relative rapidity 9 ab = iu c ah after dividing the S'-matrix by the residue 

Afabcf- 

The kinematics goes as follows. Using equ.(|25D, we can write 9 ab = iu c ab = 
i(2% — u\ a — u^) = i(v% a + «5c). This suggests to apportion rapidities as : 

e a = e + tul a} e b = e-iul C i 

where 9 is the center of mass rapidity of the afe-system. Hence we get 

Oad = 0-0d + 1*4 = + v4 a 

9 bd = 9-9 d - iut- c = + %lf h - c . 

Thus the simplest equation compatible with the above requirement is the 
bootstrap equation 

s cd (6) = s^e + iu\c)S bd (e - itpj . (27) 

It is this equation we will use as a consistency requirement to construct 
S-matrices in the next section. 

The bootstrap equation generates a powerful constraint on the possi- 
ble conservation laws 0. Namely, if we continue the two-particle state 
| A a {9 a )A b {9 b ) > to 9 ab = iu a b h will be dominated by | A c >, i.e. 

lime | A a {9 + - e/2)A b (6 - + e/2) >= f ah \ A c {9) > . (28) 

e— >U 

This has to be compatible with our conservation laws [S, P s ] =0. If we set 
6 a = 9 + iu b ali — e/2 and 6 b = 9 — vu^ + e/2, then the equation 

P s | A a {9 a )A b (9 b ) > m = ( 7s a e^ + 7s V e *) | A a {9 a )A b {9 b ) > m , (29) 

when continued analytically to the pole, yields 

P s S\A a (9 a )Ab(0 b ) >»= P s fj f\A c {9) >= l^^ s e se \A c {9) >= 

SP s \A a (9 a )A b (9 b ) > m = S^e 8 ^ + 1 b s e s6 »}\A a (9 a )A b (9 b ) > m = 

Therefore we get 

7 s e « c + 7 s e ^ = %, (31) 

valid for all f abc ^ . If all angles u c ab are known, this equation can in general 
be solved for the 7* = K a (m a ) s , . . . only at certain values of s, determining 
in this way the possible conservation laws. 
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3 Factorized 5-matrices 



Let us use the prescription of the preceeding section and construct some 
interesting S'-matrices0. Of course we still have to feed in information char- 
acterizing some particular model, such as some minimum particle content, 
symmetry properties, simplicity etc. In sections § to ^ we will link up with 
the underlying field theory. 

3.1 The ^-matrix of the Lee- Yang edge singularity 

Consider the simplest possible model, where we have only one self-conjugate 
particle A\ = Aj, which according to the bootstrap principle can be consid- 
ered as an A\A\ bound-state of itself. This means that fm 7^ and this 
will be referred to as A\ having the 3 'property for obvious reasons. Thus 
setting a = o = c= lin equ. (|3lD and using u\ x = ^ from equ. (|25"D , we get 

e -W3 + e +m,/3 = L (32) 

Rewriting this as cos(7rs/3) = 1/2, we see that this equation is satisfied only 
for s = 6n ± 1, i.e. s = 1, 5, 7, 11, 13, 17, 19, 23, 25, . . .. 

Since we have a pole at 9 = 2m/ 3 the minimal solution to the equ.(|27]) 

is 

S n {0) = F*(9). (33) 

The residue of this amplitude at 9 = 2m/ 3 equals +2i\/3. Thus, although 
the S— matrix satisfies the correct unitarity equation, the theory as a whole 
is not unitary, since the residue has got the wrong sign. The resolution of 
this problem goes noting that for the unitarity equation SS^ = 1 to hold it 
is sufficient, that the in- and out-states form a complete set. In fact from 
S\in >= \ out > and < in\S^ =< out\ we get < in\S^S\in >=< out\out >= 
1. From the completeness of the in-states we conclude that SS^ = 1 and 
similarly for S^S= /. The hamiltonian doesn't have to be hermitian 0. 

We conclude, that a non-trivial unitary theory with a factorizing S— matrix 
has to contain at least two particles. 

This model is only the first one of a whole series of non-unitary modelsQ], 
as we will show after discussing the Z(iV)-models. 

As we will see, the S-matrix equ. (p3|) belongs to the non-unitary, con- 
formal field theory describing the Lee- Yang edge singularity 0, deformed by 
the only relevant operator it contains. 
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3.2 The Ising model in a magnetic field 

Since one particle is not sufficient to produce a unitary theory, let us go one 
step further and consider an extra particle A 2 . As above let us assume that 
A\ has the 3 property and also couples to the A 2 particle, i.e. assume that 
/in, /112 and /221 do not vanish. We will furthermore assume all particles to 
be self-conjugate. 

Using now equ.(|31~D with a = b = 1; c = 2 we get 



k Xe~"" 12 + K 1 s m s 1 e +lsu ^ = K 2 s m s 2 . (34) 
Repeating the same for a = b = 2; c = 1 results in 



where 



x{ + x^= R ^, x° 2 + xr= — 1, (35) 



Xl = e m iz ; x 2 = e m *K (36) 
Eliminating the r.h.s. of these two equations, we finally obtain 

(x s 1 + xr s )^ + x 2 - s ) = i. (37) 

This is a very much overdetermined system of equations, if we have non-trivial 
conservation laws with s > 1. But, as we mentioned, miracles do happen and 
equs. (|37|) do have a solution for s belonging to a subset of integers. For s — 1, 
we get : 

X\ = exp(7ii/5), x 2 = exp(27rz/5). (38) 

This solution can be obtained making the ansatz x 2 = {x\) r . r = 1 gives 
mi = m 2 , which we don't want, whereas the above solution corresponds to 
r = 2. For s = 1 we therefore obtain the golden mass ratio 

^ = 2oos^ = ^±i = 1.6180339.... (39) 
mi 5 2 y J 



With this information, we can write equ. (|37|) as 



7T 2lY S7i /2S7T 

cos(— ) cos(— ) = cos(— ) cos(^^). (40) 
5 5 5 5 
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This equation is satisfied only for s ^ 0(mod5). Recalling that A\ has the 
4> 3 property, the allowed values for s are now 

s = 1, 7, 11, 13, 17, 19, 23, 29, 31, ... . (41) 



We may now start to use the bootstrap equation (|27|) to crank out the possible 
S- matrices of this model. Start with Su. It is supposed to have poles at 
6 = and 9 = with positive residues corresponding to the particles A\ 
and A2 in the direct channel and corresponding poles with negative residues 
in the crossed channel. Using our functions F a (9) (which we abbreviate as 
[a] ) as building blocks, we choose the following O'th order trial for Su : 

S[i } = [~] [§]. (42) 

S11 has to satisfy the bootstrap equation (G7p, whose r.h.s. becomes : 

[i] HI [i] (43) 

Here we used equ. (pO|) to maneuver the displaced ^-dependence into the in- 
dices. Due to the presence of [1/15] and [11/15] , this cannot be equal to the 
l.h.s. of the bootstrap equation. Therefore we modify S[°i to what will be 
our final version: 

SnW = [|] [|] [1]. (44, 

It is now easily seen, that the bootstrap equation is satisfied at the expense 
of having introduced new poles into Su, which depending on the sign of the 
residues, will correspond to new particles. As a matter of fact the new pole 
at 9 — — tti/15 represents a new particle A3. 

Note that up to now we have from A±A\ — > A±, A\A\ A2, A\A% — > A% 
the following results: 

1 2tt 2 2n 3 7r 



47T o 27T 



J2 



5 ' Z1 5 



< = IT- ( 45 ) 



Thus from equ. (|23f) we obtain 

mi V307 



"' :! 2 cos ( ^-"j = 1.9890437 .... (46) 
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Now we have to continue turning the crank, hoping that eventually no new 
particles will be necessary in order to satisfy the bootstrap equation equ.(|27|). 
In the meanwhile, we use it with a = b = d = l,c = 2 : 

Si 2 (0) = S n (6 - m/5)S n (6 + tti/5). (47) 

This yields 

s ^ = £ I' & & (48) 

We have a new pole with positive residue at 9 = iuf 2 — 4tu/15, corresponding 
to a new particle A4 with mass 



^ = 4cos(f) cos (|) = 2.4048671.... (49) 



m 4 t ( n\ (Itx 
mi 

The amplitude S22 can now be obtained from equ.(|2T|) with a = b = 1, c = 
d = 2 : 

S 22 {9) = S 12 (6 - m/5)S 12 (9 + tti/5). (50) 

This yields 

s 23 {0) = s n (e)s vs {0). (51) 

S 22 {9) exhibits two positive-residue poles at 6 1 = = 4^/14 and 9 = 
m\ 2 = 7u/15 , representing new particles A§ and with masses: 



mi 
me 
mi 



7T\ 




/2tt 


5y 


cos 


vl5 






f 7T 


5y 


COS 


v30 



2.957295 . . . 

3.2183404.... (52) 



We notice, that double poles occurring in, for example S 22 (Q), correspond 
to anomalous thresholds 0. 

Let us go on and construct S13 {9) = Ssi(9), putting a = b = d = 1 and 
c = 3 in equ.(j31~l). We need u\ 3 , which we lift from equ. (|25|) , using u\ x = tt/15 
and u c ab = ul a . We obtain 

S 13 (9) = S u (9 + m/30)S n (0 - ?tt/30), (53) 

which yields 

SM = [-} [-] [-] [-] [-] [-] 
13 




In a similar way we get : 

= [i] [^] [A] ([1] [|] [I]) 2 

S 33 (0) = [~] [^] [^] S 11 (0)Sv{0) (55) 

^(f) = [^] [^] [^] 2 s 12 (B)s 22 (e). 

We encounter two new particles : A 7 as a pole at 9 = 2m/ 15 in 633(0) and 
A 8 as a pole at = ni/15 in S44 (6 1 ) , whose masses are : 

772,7 
mi 

7778 
777! 

At this point we may rest, since to our deep satisfaction we encounter no 
new poles and the bootstrap program closes with 8 of them. As we shall see 
these S-matrices describe the massive field theory, obtained by perturbing the 
critical T = T c zero-field Ising model by a magnetic field. We also notice that 
an underlying structure pertaining to the exceptional Lie algebra E$ raises it's 
delightful countenance : the conservation laws equ. fl4i"]) are labeled by spins, 
which are exactly the exponents of the Lie algebra E$, repeated modulo the 
Coxeter number of E s 0. This is not completely unexpected, since the 
conformal field theory describing the critical Ising model can be obtained via 
the coset construction (E s )i®( y E s )i/(E 8 )2, where the subscript denotes 
the level of the Kac-Moody algebra. 

3.3 The Z(A0-models 

This set of interesting models are the simplest generalization of the Ising 
model, which corresponds to iV = 2. Their symmetry is actually Z(2) x Z(N), 
where the extra Z(2) factor stands for charge conjugation. They exhibit a 
very rich and interesting phase diagram^, which includes special multicriti- 
cal points, where exact solutions can be obtained, even off criticality. These 
points then become relevant for the S-matrix game. 



7T\ 






1 C ° S (36 


5y 




7T\ 






1 C ° S (1 


5y 





(56) 
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On a lattice these models are defined by spins Oi living on each site i, 
which satisfy ((Ti) N = 1 or equivalently erj = (ai) N ~ l . In the field-theoretic 
context this is translated into the property, that anti-particles are bound- 
states of N — 1 particles 0. This property requires the reflection amplitude 
to vanish, since otherwise this would correspond to non- vanishing production 
in the crossed channel. 

For the scattering of fundamental particle, we have then only two ampli- 
tudes 

< P2,Pi\S\pi,p 2 >= S u (9 12 ) (57) 

< P 2 ,Pi\ s \Pi,P2 >= S lT (6 12 ). 
Unitarity and crossing imply 

u{0)u{-6) = t(6)t(-6) = 1 (58) 

u(6) =t{in-6). 

We now make a minimality assumption and introduce a pole, corresponding 
to a two-particle bound-state at 6i 2 = iu\ x = 2m /N in u(9), the following 
mass spectrum is generated [[]]: 

sin(7ra/iV) . , 

m a = m .\ / , a = l,...,iV- 1. 59 

sm(7r/iv) 

Here we have m^-i = m or more generally, in agreement with charge- 
conjugation invariance, rriN-a = ma- 
in the present model, as opposed to the Ising case with h ^ 0, T = T c , 
we know the mass spectrum and therefore don't have to go through all the 
motions to find u c ab . We immediately obtain the minimal 5-matrix of the 
'fundamental' particle 0]: 

s 11 (e) = f>(0). (60) 



For N = 3 just apply the equ. (|27|) . For larger iV one has to fuse N — 1 
particles to check, that anti-particles are bound-states of iV — 1 particles [0], 
which means that the following identity must hold : 

\{u(6 + niit/N) =t{6) =u{m -6), (61) 
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where 

±1, ±3, . . . , ±(N - 2)forNodd 
0, ±2, ±4, . . . , ±(N - 2)forNeven. 



n 



For N = 2 this gives Su(6) = — 1, as it should for a free bosonic theory. 

The complete two-particle S'-matrix is obtained via the Bootstrap Prin- 
ciple as 



S ab (9) =iW0) 



i 2 

min(a,b) — 1 

1 J f \a-b\+2k (8) 

k = l 



f^{e), (62) 



where a,b = 1,2, . . . , N — 1. The simple poles of S a b(0), which would vio- 
late the Z(AQ-symmetry, if interpreted as particles, are doomed to exhibit 
residues with the wrong sign. 

We can now check for which values of s we encounter conservation laws. 
The mass spectrum together with equ.(^) yields 



u 



Ti(a + b) 



ab N 



c = a + b(modN). (63) 



Using this information in equ. (|31~D, we arrive at - recall a = N — a ! - 

^a^snb/N + 7 y»# = 7 *f6_ ( 64 ) 

Since the 7^ are real, the imaginary part of the l.h.s. has to vanish, yielding 
the following condition: 

7 S 6 sin(s7r&/iV) = 7s a sin(sirb/N) (65) 

This therefore better be an a, 6-independent const ant, which we normalize to 
unity and get: 

7 s a = sin(svra/iV). (66) 

The real part gives: 

7s a cos(svr&/iV) + 7s b cos(svra/iV) = 7s a+b sin[vr(a + b)/N] (67) 

and this is identically satisfied for 7° from equ.(|65D. All this is of course 
consistent only for sin(7rsa/A^) 7^ 0, or 

s ^ 0(modN) (68) 
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and we get an infinite number of conservation laws for all these values of s. 
Let us expose, that this set reveals an underlying structure corresponding to 
the group A N _ X = SU(N). 

The relevant properties of simple Lie algebras can be encoded in Dynkin 
diagrams 0. For A N _i or SU(N) this looks like 



4 



N-l 



The incidence matrix I ao of this graph is the (N — 1) * (N — 1) matrix 
I a b = number of lines between nodes a and b 

/ 1 ... \ 
1 1 ... 
1 ... 



The eigenvalues of A are 



\ 1 J 



S7T sin^f 



(69) 



(70) 



for s = 1,2, . . . , N — 1. The corresponding normalized eigenvectors are 



p(o) 



2 sail 
N Shl 'N 



(71) 



The matrix / has non-negative entries and one defines the Perron- Frobenius 
vector as the unique eigenvector, all of whose components can be chosen to 
be positive. The corresponding eigenvalue is not smaller in magnitude than 
any other eigenvalue||. In our case it is the vector Fi. 

At this point, we take a rest and compare our Z(N) mass formula equ.([59]) 
and our result for 7" equ. (|66|) with the preceeding algebraic constructs. We 
immediately notice: 

i) = sin ^ and therefore the components of the Perron- Frobenius 
vector give - up to normalization - the mass spectrum. 

ii) r^ a ) = 7° or F s = 7 S , where 7, = 7] , . . . , 7 S Y_1 . The eigenvalues of the 
conserved charges P s are thence also given by group theory. 
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If we replace in the above formulas the number N by the Coxeter number 
h A of the group A, then these two results are valid for all the simple-laced 
groups A, V, £. £ 8 corresponds to the Ising model in a magnetic field and all 
the other cases have also been identified 0, [|. 

We now want to reap some rewards by extending our solutions a bit. 
All the S'-matrices obtained in this section are mimimal in the sense, that 
they satisfy all the imposed constraints with a minimum number of poles. 
However there are models with the same symmetries, but which contain e.g. 
a free parameter. A very prominent set are the Toda field theories ||, which 
are exactly integrable and contain a coupling parameter j3, on which the S- 
matrices must depend. As first shown in ref.0, one can find solutions of 
the relevant equations, which are products of two factors : one equals the 
minimum S^^-matrix and the other Z a b(9) contains the (3 dependence. The 
mass-spectrum does not depend on (3 and the theory becomes a free one, as 
(3 — > 0. Therefore S a b(9)Z ab (9) goes to unity as (3 — > 0. If we build Z ah (Q) 
out of f a {6) with a < 0, we introduce no new poles in the physical sheet and 
are able to cancel all poles of S a b{9) as (3 — > 0. For the Z(N) Toda models 
the Zn(#)-factor for example is 

Zn(e) = Me)f_^ +fi (B). (72) 

Factorized S'-matrices describing Toda field theories with structure corre- 
sponding to all the simple Lie groups have been found recently || [§. 



3.4 A non-unitary series A'(2N) 



Finally we may generalize the non-unitary Lee- Yang edge singularity to a 
whole series of non-unitary models [[]]. This is achieved, replacing the Z(N)- 
amplitudes f a (9) by F a {6) in equ.fl62"D, such that they reduce to the Lee- Yang 
case for N — 1. In this way, we certainly obtain minimal amplitudes for a 
set of iV self-conjugate particles: 



S a b{9) = F\ a -b\ (9) 



min(a,b) — l 

n 

k=l 



a-f)|+2fc 



(0) 



with 1 < a, b < N. The mass-spectrum is given by 

G 



sin 



sin 



2JV-1 



2JV-1 



1,2, 



F(a + b) (9), 



N - 1. 



(73) 



(74) 
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As opposed to the Z(N)-models, here the absence Z(iV)-invariance places 
no restrictions on simple poles to interpreted as particles. The whole se- 
ries therefore violates unitarity i. As we will see, these S'-matrices belong 
to a series of non-unitary conformal theories with negative central charge, 
perturbed by a relevant operator. 

The investigation of possible conserved charges P s runs parallel to our 
Z(N) discussion, so that we don't repeat it here. One fmds|J, not surpris- 
ingly for us at this stage, that an infinite number of P s exists for 

s= odd ^ 0(mod(2N + 1)) . (75) 

Here also exist non-minimal solutions and their fundamental Z-factor is 
given by0 

Z ll {6) = F- P {e)F__ i _(6). (76) 

These non-mimimal models are now unitary, due to the Z-factors, which 
change the sign of the relevant residues. 



4 The Thermodynamic Bethe Ansatz 

The Bethe Ansatz (BA) has been a powerful tool in analysing 2-dimensional 
field theory and statistical mechanical models Q. It has been used to diag- 
onalize the hamiltonian of these models, to study their spectrum, finite-size 
corrections etc. Our purpose is to use the BA as a means to interpolate from 
the massive S-matrix theory to the ultra-violet, short-distance fixed point, 
in order to make contact with data from conformally invariant field theory. 

One way to use the BA is to start with the unrenormalized Hamiltonian 
and to make an Ansatz for the wave functions of 'pseudo-particles' living 
on top of a pseudo- vacuum. Their spectrum is obtained by imposing e.g. 
periodic boundary conditions on the wave functions. This yields the so- 
called BA equations. Bound states correspond to complex solution in the 
rapidity of the BA equations. 

These pseudo-particles have then to be used to fill up a Dirac sea in 
order to construct the physical vacuum. The excitation spectrum above this 
sea will correspond to the real particles of the theory, whose scattering is 
described by the ( physical ) S'-matrix. 

3 See however reference 101. 



19 



In the Thermodynamic Bethe Ansatz ( TBA ) || |(J one deals with the 
real particles from the very beginning, the rapidities therefore being real 
numbers in the physical region. Bound states are treated on equal footing, 
which is also in accord with the bootstrap principle. 

Suppose now we know the ^-matrices of all the particles in the theory. 
We may then place them in a box, let them scatter till they are in equilibrium 
at a particular temperature T. As we shall see, it is reasonable to assume, 
that these equilibrium states will be described by BA wave functions, since 
there is no production in our models. As shown Al. B. Zamolodchikov 
[[]], we may now extract finite size effects and go in particular to the zero 
mass limit. In this way we obtain e.g. the conformal anomaly ( or central 
charge ) c of the ultraviolet conformally invariant limit of our massive field 
theory, establishing a beautiful connection between factorizable S'-matrices 
and conformal invariant field theories. 

The importance of the central charge c , stems from the fact, that its 
knowledge is almost sufficient to pin down the conformal field theory we are 
talking about. As will be explained in section [5] , it can be extracted from 
finite-size effects arising in a massive Euclidean field theory defined on a 
rectangle of size L * R with , say , periodic boundary conditions in the both 
directions - i.e. on a torus. In the following, we will emphasize periodicity in 
the R-( space-) direction and talk about a 'vertical' cylinder of radius R. If 
we slice the cylinder parallel to its base, we may define the transfer matrix 
T as the partition function of the small cylinder AL * R, where we fix the 
configurations of the fields on the neighboring bases §. In order to write 
the partition function Z(L, R) in terms of the transfer matrix T along the 
L—(time- )direction, we introduce a complete set of states \n(L) > for each 
slice at position L. Due to periodic boundary conditions in the vertical 
direction, the states on the bottom and top bases are the same. Labelling 
the slices by 1, 2, . . . , 1 , we get 

Z(L, R) — i 1] El <n 1 \f\n 2 ><n 2 \f\n 3 > ... < n L _i|f|m >= 

\i=l,2,...,l n 3 J 

Y,<ni\f L \n 1 >=Trf L . (77) 

m 

4 We always assume some UV cutoff, so that we can count configurations and take the 
limit at the end. We also choose AL = 1. 
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The quantum hamiltonian of this theory is defined as T = exp(-H). Since 
the hamiltonian H is related to the time-time component of the energy- 
momentum tensor, integrated over space, we have 

H=^~ I f LL dR } (78) 

Z7T J 

where a factor 1/2-7T has been extracted for later convenience. 

If E (R) is the lowest eigenvalue of H ( the ground state energy of our 
field theory ), then for L —>■ oo it will dominate the partition function : 

Z(L, R) ~ exp(-LE (R)). (79) 

The central charge c is now obtained from the free energy per unit length 
- which equals the ground state energy E of the two-dimensional conformal 
field theory - 

F(R) = - lim - In Z{L,R) = E (R). 

L— >oo Xj 

The TBA will extract c and other interesting quantities, putting the system 
in a box at finite temperature T. 

Take as 'box' our cylinder. As L — > oo, we may look at our system 
as a Euclidean quantum field theory ( EQFT ) living on an infinite one- 
dimensional space with periodic 'time' in the indirection i . But this is 
exactly, what defines a finite temperature Gibbs state (remember k = 1!). 
Thus Z(L,R) can be regarded as a EQFT at finite temperature T = 1/R 
with free energy per unit length 

f(T) = -kT\lnZ{L,R) = (80) 

Therefore, comparing equs. (|79|) and (|S0D, we finally get the relation 

E (R) = Rf(R). (81) 

This establishes the desired link between a quantity pertaining to a EQFT 
and a thermodynamic one accessible to the TBA using S'-matrix elements as 
input. 



3 The period is R, the factor 2ir being absorbed in a rescaling of the energy. 
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As we will show in section [| equ.( |156| ), the expansion of the free energy 
( of the UV limiting theory ) in powers of 1/R is: 



F(R) = f R 



TTC 



+ 0(1/R). 



(82) 



QR 



Here /o is the non-universal bulk free energy per unit area of the infinite 
plane, which is chosen to vanish. Due to periodic boundary conditions no 
term proportional to R appears, c is related to c by c = c — 12A , where 
A is the lowest critical dimension of the conformal field theory. For unitary 
theories Ao = and c = c. In order to reach c , we still have to take an UV 
limit to get from our massive S"-matrix theory to the conformally invariant 
fixed point. From equ.(|7^) we see that E (R) has the dimension length -1 
and scaling arguments then imply 



where r = The correlation length £ is related to the mass of the 

lightest particle ( in the thermodynamic limit, where Lorentz invariance holds 
and where the momentum P = 2irn/R, n = integer becomes a continuous 
variable) as £ = TOq 1 , so that r = Rm . In the UV limit, where r « 1, c(r) 
becomes the central charge c = c(0). 

Let us now set up the TBA equations, considering a EQFT at temperature 
T = 1/R on a periodic one-dimensional space of length L. For simplicity, let 
us start considering a system of N identical particles of mass m at equilibrium 
at temperature T = 1/R. Suppose we take particle a on a round-trip along 
the circle of length L. On its way it will scatter against all other particles 
it encounters. Let us look at this process in some detail 0. If our system 
contains N particles at positions Xj, there are N\ regions in configuration 
space, where all particles are well separated - \xi — Xj\ » £ - and their 
mutual interaction can be neglected i. In each of these regions, we can 
describe the system by an iV-particle wave function. Due to a scattering 
process - made up of two-particle scatterings - , the system will pass from 
one region to another. The wave function of the final region is now obtained 
by multiplying the wave function of the initial region by the relevant ^-matrix 
( exp(— 15{6) ). On it's trip 'around the world' all the phase shifts will add 

6 Typically it will decay as exp(—x/^). 




(83) 
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up and due to periodic boundary conditions the final wave function will be 
identical to the one we started with. This has to be truefor anyone of the N 
particles. Thus we get : 

or 1 

N 

jLmsinhfr) -Q _ = 1,^=1,2,..., iV, (84) 

or taking the logarithm 

N 

Lm sinh Q i + ^ 5(6i - 6j) = 2vr J i , i = 1, 2, . . . , N, (85) 

where Ji are integers to be determined by equilibrium conditions. The above 
Bethe Ansatz equations ( BAE ) reduce to the well known equations in the 
free case o~ij(8) = 0, where they too determine the possible values of the 
momenta Pi = rrii sinh(#j) = 2nJi/L. In our non-trivial, interacting situation 
they are a set of n coupled transcendental equation for 9^. The energy and 
momentum are given by 

N N 

E = ^^m cosh 0j, P = ^2,m sinh 9j, . (86) 

i=i j=i 

Notice that factorization of the 5*-matrix once more has given us an un- 
expected surprise : how to extract off-shell information from on-shell data 

In the following we will assume our Bethe wave functions to obey an 
exclusion principle, i.e. they are of fermionic type [0, 0. All our ^-matrices 
are compatible with this condition. This seems to be a necessary requirement 
for the construction of a physical vacuum in the usual Bethe- Ansatz [[]], 
although it is not necessary for the implementation of the TBA 0. 

The wave function of identical bosons or fermions has in his case to be 
anti-symmetric under their interchange. This implies the following constraint 

7 If the S'-matrix is not diagonal, the wave functions will have internal indices, on which 
the S'-matrix acts. In this case the following equations turn into matrix equations. 

8 In another step to implement this program, one obtains the form-factors from the 
on-shell information 101. 
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on the ^-matrix. From the unitarity condition equ.([ll|), we have S' aa (0) = 
±1. If S aa (0) = — 1, the wave function is anti-symmetric under interchange. 
Therefore, in order to enforce an exclusion principle, we require S aa (0) = — 1 
for bosons. For fermion the situation is reversed and we require S aa (0) = +1 
for fermions. 



4.1 The TBA equations in the thermodynamic limit 

To obtain information about the spectrum of our theory, we have to inves- 
tigate the roots of the BAE, which is a complicated system of transcenden- 
tal equations. But this system simplifies in the thermodynamic limit, where 
L — > oo together with the number of particles, which also grows ~ L. The dif- 
ference between adjacent solutions of equ. (|85|) goes to zero as 9— 9j ~ \ jm§L. 
It is then reasonable to define a density J(6) giving the number of real roots 
in the rapidity interval 9 and 9 + A9. 

Suppose now that {9 1 , 9 2 , . . . , 9 n } is a self-consistent i solution of equ. (|85|) . 
Consider then the function 

mL N 
J{9) = — sinh0 + Y j 8{9-9 j ), (87) 

2,71 

3 

which - as will be checked later on - is a monotonically increasing function of 
9. Whenever J(9) passes through one of the integers J,, the corresponding 9 
equals Q{. However there may be integers J (9) for which the corresponding 9 
is not in the set 9i,9%, ■ ■ ■ ,9 n . Such a 9 will be called a hole. Thus a natural 
definition of the density of roots and holes is 

A linear integral equation for J (9) can now be obtained in the limit 
L — > oo, by first differentiating equ.(|87|) 



1 d J{6) - p{9) = ^ cosh£ + £ V {9 - (89) 



L d9 rw 2tt 



9 By this we mean that the 6j, besides obeying the BAE, also satisfy the equilibrium 
conditions at temperature T to be set up below. 
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where <p(9) is defined as 

<p(e) = ±6{ff) = -t±hxS{ff). (90) 
We now make the replacement ( valid for sufficiently large L ) 

EM)- J de P {e)f{e)-Y J f{e j ) (91) 

d9 Pr (9)f(9), 

where 61, . . . , 9\ are the hole rapidities and we introduced the density of roots 
p r {9). Clearly the hole density is: Ph = P — Pr- 
This transforms equ.(^) into 

Tfl 

p(B) = — cosh 6 + if* p r {6), (92) 

27T 

where we have introduced the convolution 

(f*p)(9)= [ +00 ^-f(9-9')p(9'). (93) 

J -00 Z7T 

In general we may have more than one type of particles. Therefore in- 
troduce densities for each type, labeled by a — 1, ...,n : p( a '(9), etc. The 
integral equation for p^ = p^ + pfi' becomes then : 

p ("\6) = ^cosh6 + jh Vab * p?\6), o = l,. ..,n. (94) 

2n 6=1 

As advertised, the very hard problem of solving N coupled transcenden- 
tal equations has been converted to the solution of a simple linear integral 
equation for the density pf 1 '. However we still need a relation between p)?' 
and pjfi. 

In order to obtain the necessary additional information, we impose the 
condition of thermodynamic equilibrium. Thus we will calculate the free 
energy per unit length f = e — Ts, whose minimization will provide what we 
want. 
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It is easy to compute the entropy of the distribution of particles and holes 
on the circle L. We have a large number of roots Lp^ AO and holes p^ A9 
in the interval AO. The number of ways to distribute them in this interval - 
according to our exclusion principle - is 



l(pP + pi a) )(e) Ae 



(95) 



[L P f\e) Ae]\[L P t\e) ao}\ 

Taking the logarithm, we obtain the entropy per unit length: 

n 

S \piPr] = E S a[p,Pr] 
a=l 

p ~\~ oo 

= E / M [(P { r a) + HP [ r a) + Pt ] ) ~ Pl a) M a) - pi a) lnpi a) ]. (96) 
The energy per unit length is obviously given by 

n r+oo 

e[P,Pr] = E / d6pi a \6)m a cosh 6. (97) 



a=l ' 



Let us introduce the function e a {6) by 
and also 

L a {9) =ln(l + e- £a(e) ). (99) 
Using this the entropy may easily be expressed as 

n /"+00 

s[p r ,e a } = Y: / de P ^\e)[tM + ^ + e^ e) )Hi + e-^)]. (100) 

a=l J -°° 

The equilibrium distribution at temperature T is now given by minimizing 
the free energy f = e — Ts with respect to p*f\ subject to the periodic 
boundary condition equ. (|94|) relating p^ to p[ a \ 
Thus using 

5s 5s 5s 5p 

5p r 5p r 5p 5p r 

= In^ + ^ln^- (101) 

Pr 5p r P~ Pr 
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and 

coming from the periodic boundary condition, we get finally get the following 
extremum condition: 

m n 

-^coshO = e a (6)+^2fab * L h {0). (103) 

1 6=1 

Note that this condition depends only on e a (9). The solution of this equation 
for e a (9) will always be real. From equ.(^) we see that this implies p^ a \6) > 
0, since p < f\9) is positive by definition. This in turn results - see equ.(^) - 
in j( a \9) being a monotonically increasing function of 9, as announced. 

The equilibrium free energy can now be obtained by eliminating p^, using 
the periodic boundary condition equ.(p4p and p^j p r = e e ^ in the expression 
for s[p r ,ph\. Using the extremum condition equ. (|103|) then yields the free 
energy : 

n r+°° d9 

f{9) = / —m a cosh9L a (9). (104) 

a=1 J-oc 2tt 

Usually numerical work will now be envolved in order to compute this equilib- 
rium free energy. But a number like the conformal anomaly can be extracted 
analytically. 

To accomplish this we have to go to the UV limit, where the masses may 
be neglected : r = m/T — > 0. Let us first take the limit of the extremum 
condition equ.( 



m, 



a rcosh9 = e a (9) + J2<Pab * L h {9), (105) 



6=1 



where we have introduced dimensionless mass ratios m a = m a /m = m a £. 
As r — > 0, 9 has to go to oo for the l.h.s. to give a finite contribution, in 
which case it behaves as 

rh a r cosh 9 ~ rh a -e = m a exp(# — In -). 

2 r 

We therefore make the shift 9^9 + ln(2/r) to get the following equation: 

n 

e 9 = e a W + E(^ * Lb){9), (106) 

6=1 



27 



whose solution provides the r -independent functions 



e o (0) = e(0 + ln(2/r)) (107) 

L b (9) = L b (9 + \n(2/r)). 

Solving numerically the extremum condition for e a (9), one realizes that this 
function equals the constant e a in the range — ln(2/r) << 9 << +ln(2/r) 
and eventually grows exponentially, so that e a {9) = e a , except at oo, for 
r->0 |. Similarly L b (9) interpolates between e a at 9 = and - through a 
double exponential decay - for 9^0. 

If we now use equ. (|83|) and take the limit r — > in equ.( |105| ), we get a 
compact formula for c(0): 



C (°) = ^E^/ d9L a (9)e\ (108) 

^ o=l J -°° 



which is very convenient for a numerical solution. 

A perhaps more explicit formula for c(0) may also be obtained in terms 
of Roger's dilogarthmic function [[]|. It is more convenient to use the entropy 
equ.( |100|) and to get the c(0) from 

S = ^P-TL + 0(T 2 ), (109) 

which is equivalent to equ.(|8"2"D 0. 

Again the entropy s[p r , e a ] vanishes in the limit r — > 0, as long as the limits 
of the integral over 9 are finite. Now check the behavior of the integrand. On 
one hand taking the derivative with respect to 9 of the extremum condition 
equ.( |103|) we get 

de a (9) m a m a m 

~ — sinhfe* ~ — e l l sqn9, 6* — > oo. 110) 

dQ T 2T K ' 

On the other hand, since <p{9) = + 0(e~ 8 ) as 9 — > oo, we get for the 
asymptotic behavior of from the BAE equ. flP]) : 

Air 



10 We have f(T) =e-Ts = T 2 (tt£/6). But s = -df/dT = T(7r5/6) + 0(T 2 ) and 
equ. dlO|) follows. 
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and consequently 



pl a) * 



: ,\e\ 



iTi 1 + e Ea ^ 



T de a (6) 



1 + e 



(sgnd). 



2tt dd 

Substituting this into equ.( |100|) for the entropy, we get 

T A f+~ .MB),,,.,, , . m , £„(9) 



1 + c^W) 



(112) 



|. (113) 



Here we change variables 6 — * e a and use ip ab (9) = <p a b(—@), which follows 
from the unitarity of the S-matrix , to show that e a (6) is even in 9. Finally 
changing variables again from e a to /(e a ) — 1/(1 + ex P( e a)), we get : 



lim s\p r ] 



T » /-/M°o)] , r lny ln(l-y). 
2tt ^J/ie.to)] 1-y y 



(114) 



The upper limit of the integral equals 0, since e(8) — > oo there. The lower 
limit may be obtained from the extremum condition equ. (|103| ) , setting there 
6 -> and r -> 0: 

= e«(0) + £ / — ^(0 - 6')L b (6'). (115) 

However ip ab (9) decreases exponentially off the origin and we get an equation 
for e a = e a (0) 

(116) 



e e « = J|(l + e ~ e ») Nab 

6=1 



a = 1, . . . , n, 



where the symmetric matrix N ab is: 

r+oo d6 1 
N a b = ~ TT^abiO) = - — [5 ab ( + (X)) - S ab (-00)}. 

J-oo ATX Z7T 

Introducing Rogers' dilogarithmic function L(x) as 



(117) 



L(x) 



dy 



lny ln(l - y) 



1-?/ 



(118) 
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we may finally express c, using the finite-size scaling formula equ.( |109| ), as 

n 

c = ^c a (e a ), where 

a=l 

*) = ^ L (TT^) = >( e i < 119) 

Each particle species contributes with c a (e a ) to the total central charge. The 
function c a (e a ) is strictly monotonically decreasing, approaching as e — > oo 
and c Q (0) = 1. In order to obtain c the non- linear equation ( |116|) has to be 
solved. 



4.2 Computation of the central charge 

Let us try things out in some simple cases, although there exist solutions for 
whole sets of models 0. 

1) Ising model ( h = 0, T ^ T c ) 

The Ising model ( or Z(2) model ) is a free fermion theory with <^ a & = 
const and N a b = 0. We may start directly from equ.( |104|) for the free energy. 
Thus we have to compute 

E (r) = Rf(R) = -^c = 



/+oo rjf) 
—r cosh#ln(l + e- e{e) ) 
-oo 2tc 



1 r +oc de 
~R 

in the limit r — » 0. With e(9) = rcosh^ ~ re' 9 '/2, we get 

/ + 0O rlfj fOD rlfj POO 

— cosh#ln(l+e- e(e) ) ~2 / — e{9)' ln(l+ e - e(e) ) ~ / de ln(l+ e - £ ) 
-oo 2tt Jo 2ti Jo 

With this result we obtain, as expected c = 1/2. As a matter of fact, we 
may obtain an expansion of RE (R) in powers r and the result compares 
well with numerical finite-size calculations [[J. 

2) Lee- Yang edge singularity [|| 

In section |[ we found the 5-matrix 

Sn(e) = F 2/3 (e). (120) 
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From the properties equ.(pOD, we have 



F a {6) = f a (9)f a (m -6) = f a {6){-l)h{e)h- a {6). 



Since all our S'-matrices are products of the / a (0)'s, their contributions 
to the phase shifts and the matrix N a b of equ.( |117| ) has the structure N a b = 
Ej N[f aj }. For N[f a ] we get N[f a ] = (1 - \a\)sgn(a) for 1 < a < 1 and 
sgn(0) = 0. 

Consequently we get N[F 2 / 3 ] = 1 and equ. (|116|) becomes 



1 + e 



121) 



whose solution is 



ei = log[(V5 + l)/2] (122) 

and c can be extracted form equ. (|119[) to give 5 = 2/5. If we take the value 
of Ao from conformal field theory: Ao = —2/5, we get for the central charge 
c = —22/5. Notice the negative value for c, as appropriate for non-unitary 
models. 

3) Z(N)-models 

Now let us look at Z(N)- (also called Ajy-i) models in some detail. The 
complete two-particle S'-matrix was written in equ.(E^): 



S ab {9) = {\a-b\/N} 



min(a,b)—l 



n { 



k=l 



\a-b\ +2k 
N 



{(a + b)/N}, (123) 



where a,b = 1, 2, . . . , iV — 1 and we used the shorthand {a} = f a {d). 

Let us check Z(3). We have : N n = N 22 = N[f 2/3 ] = 1/3 and A^ 12 = 
AT 21 = N[f 1/3 ] + N[f 3/3 ] = 2/3. Therefore we have to solve the equs.([TT6|): 



e £1 = x\ 



e 62 = x 2 



(l + e- ei ) 1/3 (l + e~ e2 ) 2/3 (124) 
l + e -i)2/3 (1 + e -e 2) i/3_ 

Their solution is x\ = x 2 = 2cos(7r/5). If we now use the sum-rule [[J 

2n it 2 2n 



^ ( sin 2 (vr/n + 3) 

k=2 



sm 2 (7:k/n + 3) 



L(l) 



n + 3 6 n + 3 



(125) 



31 



we obtain c from equ.( |119| ) as 

6 2 6 7r 2 4 4 

TT 2 i 5 5 



which is the expected value c(Z(3)) = 4/5. 
For general A" the solution of equs. ( |116| ) is 



air ■ (a+2)ir 
Sm lV+2 S111 N+2 

sin 2 * 



(127) 



7V+2 

and the same sum rule plugged into equ.( |119j ) yields : 

c ( Z ( N » = 2i N + 2 ] > (128) 

which is the central charge of the Z(N) parafermion models. 

The non-minimal models describing Toda field theories have, as shown in 
section [3|, a free UV limit and therefore their value of c is simply the number 
of bosonic fields, as can be explicitly verified || . 

For further applications of this scheme see e.g. ref. ||. 

3) The non-unitary series 

Following the reasoning exhibited for the Lee- Yang edge singularity and 
the Z(N) models, we obtain: 

e £a (A'(2N)) = e ea (Z(2N)) . (129) 

From the same sum-rule used in the Z(N) case we get: 

~c(A'(2N))= 1 -c(Z(N)) = ^^. (130) 

Borrowing A = — ^^n+i^ > we ^ or ^ ne cen t ra l charge: 

/ .// u 2(N-1)(6N-1) , . 

c (A'(2N)) = - 1 2 J K +1 1. (131) 

At this point we have completed the first of our connections between a 
massive ^-matrix theory and it's UV-limiting fixed point data. It gives us 
confidence, that the plausibility assumptions made in the previous section, 
are indeed correct. In order to further clarify this point, we will dig into this 
UV theory in more detail in the next section. 
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5 Lightning Overview of Conformal Invari- 



ance 

The purpose of this section is to provide a minimum of vocabulary for the 
newcomer or a reminder for the rusty practitioner, but it is not a substitute 
for the excellent review articles and books on the subject |], |(J. 

5.1 Conformal transformations 

Conformal transformations are coordinate transformations x M — > y^x), which 
preserve the angle between two arbitrary vectors at any point, although their 
length may change: 

dy,dy» = ^^dafda? = p(x) dx,dx». (132) 

A special case are dilatations, which scale all lengths up by a factor p(x) = A. 
Under mild assumptions - for D=2 a discrete spectrum of dimensions - a scale 
invariant quantum field theory is also conformally invariant [[J. Therefore 
in Statistical Mechanics the long distance behavior at a critical point, which 
is scale invariant due to a diverging correlation length, is described by a 
conformally invariant Euclidean field theory 0. 

Although for dimensions D > 2, the conformal group is finite-dimensional, 
for D = 2 it involves an infinte number of parameters, which entails powerful 
restrictions on the theory. 

Consider then a Euclidean Field theory in the 2-dimensional space of 
coordinates Xi,X2 0. We, very conveniently, introduce complex ( light-cone 
) coordinates 

z = x 1 + ix 2 ; z = x 1 — ix 2 . (133) 
Then any transformation 

z ^ z ' = f( z ) ; z^z' = J(z), (134) 

where f(z) and f(z) are differentiable functions, are conformal transforma- 
tions satisfying equ.( |132j ). We will consider z and z to be independent com- 
plex variables. However, since physical correlation functions live in the real 

11 The Minkowskian space-time would have coordinates x 1 and x 2 = it. 
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space R 2 , we eventually have to take into account that z is the complex 
conjugate of z: z = z*, when constructing physical correlation functions. 

Under the finite conformal transformation z — > z' = f(z),~z — ► z' = 
f(z) the Jacobian is 

Wi= nz)m - (135) 

Tensor fields will transform with some powers of /' and / . We will define 
conformal fields of weight (h, h) to transform as 

4>{z,z) -> <j>'{z,z) = f'(z) h J'(zfy(f(z)j(z)). (136) 

The scaling dimension is A = h + h, whereas the spin is s = h — h. We also 
want to consider f(z) and f{z) to be an arbitrary, but infinitesimal change 
of coordinates 0: 5z = e(z), Sz = e(z). The above equation then becomes 

5<j)(z,z) = (e(z)d z + he'(z) + e(z)c\+he'(z))(f)(z } z). (137) 

and fields transforming this way are called primary. All fields we will deal 
with are primary, except the energy-momentum tensor. 



5.2 The energy-momentum tensor 

The energy-momentum tensor T^ v is the generator of space-time symmetries 
and as such is symmetric and conserved: T^ v = T UfJ- , d^T^ u = 0. Due to these 
properties it is a prominent object in the study of conformal transformations 
in field theory. If we consider the infinitesimal coordinate transformation - 
not necessarily conformal ! - 

x» x» + e'*(x), (138) 
then we may define T^ u as inducing the following change in the action 111: 

5S=~ f T^e v d 2 x. (139) 

Z7T J 

12 In doing this we want to stay within the conformal plane. Thus by finite we mean 
global transformations of the form z' = ffq^ with ad — be = 1, which are the only 
one-to-one mappings of the complex plane onto itself. 

13 Let us mention a delicate point here. Infinitesimal, everywhere analytic functions do 
not exist, since they are bound to develop singularities somewhere, if they are nontrivial. 
Thus we consider a bounded domain, which excludes the singularities. 

14 For example, in the case of a free massless boson with Lagrangian density C = h [d^) 2 , 
we get with this definition 7),„ = —d^dvcj) + \8 tiu d a 4>d a (f). 
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In terms of correlation functions this is equivalent to the identity 



n 

< <fa.(Xl) ■ ■ ■ $e<t>j{Xj) ■ ■ ■ <frn(Xn) > = 

72 

" / TT d M x ) < T^(x)Mxi)---Mxn) >, (140) 

J Z7T 

where 5 e <p(x) is the variation of <f>(x) under the transformation equ.( |138| ) 
El. This is a very useful definition of T^, if we don't know the Lagrangian 
of the theory we are studying. If we are after exactly integrable models it 
is often more profitable to shift the emphazis to other types of structures 
like symmetry properties, operator product expansions, fusion rules etc.7 the 
corresponding Lagrangian being determined a posteriori. 

If we introduce the fields T = (T 11 - T 22 + 2iT 12 ) /2 and T = (T 11 - T 22 - 
2zT 12 )/2, then the conservation equations T Miy = T UfJ- , d^T^ u = become 

&zT = ~d z &/2 ; d z T = -dzQ/2; (141) 

where G = T 11 + T 22 is the trace of the energy-momentum tensor. In a 
conformally invariant theory this trace vanishes and T and T are analytic 
functions of z and z respectively: 

T = T(z), T = T(z). (142) 

The 2-dimensional problem is thus reduced to effectively two one-dimensional 
problems and it is this feature, which makes conformally invariant models 
exactly soluble. Equs. fll42"D mean that the correlation functions 



< T(z)(f> 1 (zi, Z!) ... <fr n (z n , z n ) > (143) 



15 Fast derivation ! Define correlation functions by the functional integral (or Gibbs 
average on a lattice in statistical mechanics ): 

< (pyx) . . . >= 



JV(P 



e -s[, 



In the numerator of the functional integral on the r.h.s. change the dummy variable <p(x) — > 
(p{x)+8<:(p{x). Expand everything to first order in 8 e (p(x), getting two terms: one from (p(x) 
and another from the variation of the action S[<p + S e <p] = S[(p] + J d 2 x / (2tt)T^ d )Jl e u (x) . 
The zeroth order term cancels the l.h.s. and we get the above Ward identity. 
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are single- valued analytic functions with singularities (at most poles of finite 
order) only at the points zi, z%, . . . , z n and similarly for T(z). 

To see this integrate the r.h.s. of equ.( |140| ) over a domain excluding small 
regions around the points Xj and choose for e(z),e(z) functions vanishing 
sufficiently fast at infinity to allow partial integrations. In this domain d^T^ 
will be zero and only the surface terms around the points Xj survive. Putting 
the origin at one such point, and choosing for the excluded region a small 
circle of radius p, yields the surface term / da^T^e^. With 



do a = d 2 x 9F } X ^ 5(F(x ) - p) , Fix) = x\ + x\ 
ox** 



we get 

J da^ = j dOx^. 

Transforming from rectangular to complex coordinates, we obtain for the 
r.h.s. of equ.( |140| ) the result: 

J da^e v = ^UdzT(z)e(z)+fdzT(z)e(z)Y (144) 

Thus the Ward identity for T(z) (and a similar one for T(z) ) becomes : 

n 

< <f>i{zi,zi) ■ ■■54>j{zj,z j ) ...>= 

n 

53 { e {zj)d Zj + hj^(zj)j < . . . <p n (z n ,z n ) >= 

3=1 

f dz 

f — e(z) <T(z)(j)i(z 1 ,zi) . . .(f> n (z n ,z n ) >, (145) 
where C encircles the points Zi, . . . ,z n once in the positive sense. This state- 



ment is equivalent to - as can be seen multiplying equ.( [L46| ) by e(z) and 
integrating along C - 

< T(z)(f)i(zi,Zl) . . . (j) n (z n ,Z n ) > = 

n ( h- 1 \ 

( _ J ^ 2 + 9 *j < M z U z l) • ■■4>n{z n ^ n ) > (146) 

j=l \{ z z j) z z j ) 
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and a similar equation for T(z). These relations express the transformation 
properties of primary fields under conformal transformations. 

Surprisingly the energy-momentum tensor itself is not a primary field. As 
a matter of fact, in a non-trivial theory the two-point function < T(z)T(0) > 
cannot vanish and since the dimension of T(z) is h = 2 @, scale invariance 
implies 

< T(z)T(0) >= ^. (147) 

The dimensionless constant c is called central charge or conformal anomaly 
and it's knowledge characterizes to a large extent a conformal model, c has 
been normalized such that the free scalar massless boson has c = 1. The 
analog of equ.( |146 | ) for T{z) itself now becomes 

™ = (i^r + v^tf™ + J^) 9lT{z '- ] + ■ (148) 



Equ. ( |148|) translates into the following transformation law for T(z): 



5T(z) = ^ce'"(z) + 2e'(z) + e(z)d z T(z). (149) 

It is not completely trivial to integrate this quation to obtain the trans- 
formation law for finite conformal transformations. We state only the result 

T( z ) = (f(z')yT(z>) + ^S(f,z), (150) 
where the Schwartzian derivative is 



S(f,z) 



This is a very useful result, for it permits to connect informations per- 
taining to different geometries. For arbitrary f(z), not of the form f(z) = 



"This can already be seen from equ.([146|). Or : the generator of translations by a vector 
a M is exp(za M P M ). Therefore P M = jT^da^ has dimension one and has dimension 
two. 

17 This expression is to be understood as occurring inside a correlation function. The 
ellipsis stands for terms, which are regular as z\ — > generated when T(z\) hits other 
fields of the correlation function. The reader, who is uneasy with this kind of operator 
statements should take refugee in Furlan et al.Oj. 
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(az + b)/(cz + d), the conformal plane will be mapped into a different ge- 
ometrical domain. We now broaden the concept of conformal invariance to 
include invariance under this active transformation. Let us, for example, 
conformally map the complex z-plane ( without the origin ) onto a periodic 
horizontal strip of width R in the w = u + w-plane by z = exp(2irwR) and 
analogously for z. We now postulate that objects in the z-plane go, via this 
mapping, over to the corresponding objects in the strip. Applying this to 
the energy-momentum tensor, we get 

T strip (u) = A 2 [T plane (z)z 2 - c/24]. (151) 
K 

In particular, if we now take the expectation value of T(z), since < 
Tpiane(z) > has been renormalized to zero, we find 

c 2vr 2 . , 

< T stnp (u) >= -—— . (152) 

This formula allows us to measure the central charge c exploring finite size 
effects in the statistical mechanical version of our Euclidean quantum field 
theory, where c is obtained via the free energy. 

In fact the variation of the free energy is given by a formula analogous to 
equ.( |140|) with no fields on the l.h.s. : 



f d 2 u 

5lnZ = - —— < T^(m,m) > d^e,,, (153) 

JV Z7T 

where the integration domain T> is the infinite strip. Choose a transversal 
dilatation by 5r] of the strip : e\ = 0, e2 = u\8r]. We get T^ y d^e u = (T(u) + 
T(u))Srj. The accessible statistical mechanical observable is the free energy 
per unit (vertical) length with suitable boundary conditions in the vertical 
direction 

F(R) = - lim ±-lnZ{R,M). (154) 

Af— M30 JVl 

Under the transversal dilatation F(L) changes by SF(L) = R5r] dF(R) jdR 
so that from equ. (|153|) follows, that 



dF/dR= (2n/R 2 )(c/ '12). (155) 

Using F(oo) = 0, we get 

F(R) = -c|l (156) 
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This equation holds true, if all operators have positive dimensions, as is the 
case of unitary models. As we shall see, for non-unitary models, we have at 
least one operator with negative dimension A < 0. Then equ. ( |15(j| ) has to 
be modified to c — > c = c — 12 A . This equation is then used to extract 
the value of c by numerical finite-size studies or, as we have done, using the 
Thermodynamic Bethe Ansatz. 



5.3 The representation space of the Virasoro algebra 

In order to construct a Hilbert space, in which the conformal fields may act 
as operators, we have to go from correlation functions to an operator for- 
malism. Remember that the charges, which via equal time commutators, 
genenerate the variations of the fields under some transformations, are writ- 
ten as integrals over space-like surfaces. Looking at equ.( |144| ), we see that 
this integral is in our case § dzT(z). This is made explicit adopting the radial 
quantization in the complex plane. One considers a Minkowski space a , a 1 
with a periodic space coordinate a 1 . This amounts to quantize the theory on 
a cylinder of radius R, which then acts as an infra-red cutoff for the mass- 
less fields. These separate into left- and right-moving fields living on the 
light-cones cr° ± a 1 . With the map z = e a ° +ia ' , which we have already used, 
the cylinder is then mapped onto the complex z-plane and analogously for 
right-movers. The space-like 'equal time' slices are concentric circles around 
the origin, whereas the 'time' direction is radially outward, t = — oo corre- 
sponding to z = and t = +oo to the point z = oo of the complex plane. 
Time translations cr° — > cr° + r on the cylinder are mapped into dilatations 
z — > ze T ( and z — > ~ze T ) in the plane. Therefore, what we call Hamiltonian 
on the cylinder, is the generator of dilatations (= scale transformations) in 
z and z. 

Remembering that Euclidean correlation functions correspond to time- 
ordered Green-functions, we can consider operators <fi(z) defined such that 

< Mzi) ■ ■ ■ (t>n{z n ) > = < 0\1l [4>x{z X ) . . . 4> n {z n )) |0 >, (157) 

where we introduced the radial ordering operation 1Z: 

= { itlf^l \ Zl \ l M (158) 

I <P2{Z2)<Pl{Zl) \Z 2 \ > \Zi\. 
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With this understanding, we can write equ. ( |145| ) as 

5$(z,z) = — I { -I ) dwe(w)R(f(w)cf)(z,z), (159) 



2m \J\w\>\z\ J\w\<\z\ 

he difference of the two line integrals giving the integral around the point z 
This is now an operator statement. 
For the particular case of T(z) itself, we get 

Sf(z) = ^ce" + (2e'(z) + e{z)d/dz)f{z) 



e{w)R(t{z)T{w)) . (160) 



l\w\>\z\ J\z\>\w\ 

We now expand T(z) and T(z) in a Laurent series: 

+oo T A +oo ~T~ 

T{z) = E n-z) = £ (161) 

n=— oo n=— oo 

where the factor z -2 has been introduced so that the operators L n have 
dimension n t3. If this expansion, together with a similar one for e(z), is now 
inserted into equ.( |160| ), we obtain the infinite-dimensional Virasoro algebras: 

[L n , L m ) = (n - m)L n+m + —cn(n 2 - l)8 n+mfi 

[L n , L m \ = (n - m)L n+m + -^cn(n 2 - l)5 n+m , 

[L n ,L m ] = 0. (162) 

The space in which the representation of the Virasoro algebra acts is then 
constructed as follows. We want to build up a Hilbert space in which the 
primary field <ph{z,z) with conformal weight (/i, h) may act. Define the vector 

\h >= ]jm<f> h (z,z)\0 > . (163) 

z^O 



18 This is what we mean by the equal-time commutator [J dxA(x), B(y)]. For euclidean 
Green functions, the time evolution converges only for At > in exp(—ArH). Therefore 
the 7?.-operation does give sense to equ. (159). 

19 We may expand be around any point, not necessarily around the origin as done here. 
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From equ.( P-46|) it follows, that 

[L n , <j>{z)\ = I ^w n+1 T(w)(f)(z) = h(n + l)z n 4>{z) + z n+1 d(f)(z), (164) 
J 2m 

consequently [L , 0(0)] = h<f)(0) and [L n , (f)(0)] = 0, n > 0. The state \h > 
therefore satisfies 

L n \h >=L n \h >= for n>0 (165) 

L \h>=h\h>, L \h>=h\h>. 

From the commutation relations of the L n 's, we see that L n decreases the di- 
mension by n units, so that the above equations guarantee, that the spectrum 
of dimensions is bounded from below. 

Notice that L_ n , L_ n for n > create new states and therefore act like 
creation operators, whereas L n , L n are destruction operators. The space is 
therefore built up by all vectors of the form 

L_ m . . .L- nN L_ mi . . .L_ mM \h > with nj,m k > 0. (166) 

These vectors are also generated acting with descendant fields on the 
vacuum. These fields are the regular terms in equ.( |146| ) as z — > w in the 
following equation, which we now write out explicitly : 

T(z)(f)(w,w) = J2( Z ~ w) n - 2 L_ n (f)(w,w) 

n>0 

1 1 _ 

L (f) + L_ l( f) + L_ 2 (f) +(z- w)L-z4> + .... (167) 



(z — w) 2 z — w 

The descendant fields are 

r cl 7 ! 

0(-«)( w ,«y) = L„ n( f)(w,w) = f —- T (z)(f)(w,w). (168) 

J w zm yz — w ) 

Comparing with the l.h.s. of equ. (p6|) , we note that (o) = L o = h(f> and 
</>( _1 ) = = 9^00. The most important descendant field - and therefore 
not a primary field ! - is the energy-momentum tensor, which is a level two 
descendant of the identity: 

(L- 2 I)(w)=£^—T(z)I = T(w). (169) 
J 2m z — w 
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This fact allows to write down differential equations satisfied by correlation functions. 
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From equ.( |168| ) and the definition of L n as 

we see that the representation space is generated by the descendant fields as 

L_ n \h >= L_ n (0(O)|O >) = (Z_ n 0)|O >= <^(0)|0 > . (170) 
This representation is irreducible, unless there exists a null vector 

\Xh+N >= ^2 • • • Lk m \h > 

with kj — N, satisfying 

L n \Xh+N >= 0,n > L \xh+N >= (h + N)\xh+N > (171) 

for some positive integer N. 

For example for N = 2, which is the first non-trivial case, these conditions 

are 

L +1 {L^ + aL 2 _ x )(j){z,z) =0 (172) 

L +2 (L_ 2 + aL\)(j){z,z) = 

and conditions for n > 3 follow automatically from the Virasoro algebra. 
Now we move the destruction operators L +1 ,L +2 to the right till they give 
zero due to equ. (|171|) . The resulting two equations are 

-- 3 „_2M5- 8ft ) 



2(2/i + l)' 2/i + l 

We will use analogous equations for level 3 degeneracy in order to obtain a 
first nontrivial conservation law in the next chapter. 

Since the subspace generated by applying L n ,L m 's with n,m < is in- 
variant, it has to be factored out in order to get an irreducible representation. 
This factor space is called a degenerate irreducible representation space (or 
modul ) generated by the degenerate field (j)(z) and iV is it's level. 

These degenerate fields are very important, because for them the operator 
algebra 

<5>l$k = Y, C lkn$n (174) 
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This equation is written in an abbreviated notation only. 
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closes with a finite number of terms. The corresponding central charge can 
be labeled by two positive integers p, p' with no common divisor: 

c(p,p') = 1 - ^ZJH. (175) 

For each value of c(p,p') there are (p — l)(p' — l)/2 primary fields with 
dimensions 

h -h - ( r P- s P r ) 2 - (p-pO 2 ri7fi x 

(i r ,s — n>p'-r,p-8 — App') v 1 ' / 

1 <r 1, 1 < s <p-l. 

For p' — p > 2, there will exist negative weights corresponding to growing 
correlation function in non-unitary theories. For c < 1 all unitary theories 
are special cases of the above with p' — p = 1 or : 

c(m) = l- 6 , m = 3,4,... (177) 
m(m + 1) 

For each value of c(m) there are again m(m — l)/2 allowed values of h: 

h rs (m) = h m - rm+ i- s = — — —, —r (1^8) 

4m(m + 1) 

with the integers r, s satisfying 1 < r < m — 1,1 < s < m. 

We also know that the representation with highest weight h rs is degen- 
erate at level rs. Besides this the correlation functions of these fields satisfy 
linear differential equations of order rs. 

We may now clarify the shift c — > c for non-unitary models. Recalling 
equ.([78D of section |], the quantum hamiltonian H is given by: 

H= ^ [ R T uu (v)dv = f R (T(v) + T(v))dv = 

Z7T JO 271 JO 

T (L» + L )--. (179) 

As advertised, we see that time evolution on the cylinder, corresponds to 
dilatations on the plane. From equ.([7^), it follows that states with higher 
eigenvalues of L and L , contribute exponentially smaller corrections to 
Z(L, R). However in non- unitary theories some operators may have negative 
dimensions. Suppose there is one such operator with dimension Ao < 0. In 
this factor exp(— 27tAq/R has to be kept in Z(L, R), effectively shifting 

c to c = c — 12A . 
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5.4 Characters of the Virasoro algebra 

For Zamolochikov's counting argument, which permits a partial identification 
of the infinite number of conservation laws surviving perturbations breaking 
conformal invariance, we need to determine the dimensional decomposition 
of the representations of the Virasoro algebras. In the representation space 
of highest weight h the characters Xh{q) are the generating functions for the 
number of linearly independent vectors at level n, therefore having eigenval- 
ues h + n of Lq. For q, \q\ < 1 we define 

oo 

Xh (q) = q-^Tr h q L ° = q h ~ c / 2A £ d h (n) q n . (180) 

n=0 

Here dh{n) counts the degeneracy of the states in the representation at level 
n. An analogous definition holds for the right Virasoro algebra. 

If there are no null states in the representation of weight h, the states at 
level n are of the form 

k 

L- ni L_ n2 . . .L-n k \h > J2 n i =n - ( 181 ) 

i 

In this case dh(n) equals p(n) - the number of partitions of the integer n. 
Euler's generating function for these partions gives 

OO OO 1 

E * n P(n) = 11(1 - gT 1 - 1/24 „ (182) 

n=0 n=l y '/W 

and 

Xk(q) = q h+{1 - c)/2A v(q)- 1 - (183) 

Here we have introduced Dedekind's //-function, central in the study of elliptic 
functions. The c-dependent factor has introduced, so that Xrs(o) has nice 
modular transformation properties. 

Let us compute the characters for the degenerate, unitary series with 
c < 1. The allowed weights are given by equ. (|178|) and are degenerate at 



level rs. Thus in counting the states we have to subtract the null state at 
level rs and all its descendants, getting 

xrs{q) = q^ c)/2i v{q)-\q hr < s - q hr < s+rs + •••)• (184) 



44 



But the null state has weight h r ^ s + rs = h r ^ s = /i m+r . jm+ i_ s and therefore we 
have to subtract in turn it (and all its descendants) at level (m+r)(m+l — s). 
Correcting the above formula for Xrs, we get 

Xrs (q) = g( 1 - c )/ 24 r / (g)- 1 (g /l '^ - g^.-(l - q (m+r)(m+l~s) + _ ^ 

= g (1 " c)/24 r/(g)- 1 (g^ fl - q h <~--° + q h ^+r, s _...). ( 185 ) 
Repeating this process yields finally the correct expression for the character: 

oo 

Xrs(q) = g (1_c)/24 r7(g) _1 £ (^-^-). (186) 

fc=— oo 

6 Deformations of Conformal Invariant Field 
Theories 

In this section we will study relevant perturbations of Conformal Invariant 
Field Theories ( CFT ), i.e. perturbations which drive the system away from 
its UV critical unstable fixed point. An example would be the perturbation 
of a statistical mechanical system at its critical point by a magnetic field, 
the raising of the temperature off T = T c etc. The system will then flow 
away from its UV fixed point and may either end up at another critical 
conformally invariant fixed point or develop a finite correlation length, i.e. 
the theory becomes massive. 

In a CFT conservation laws are trivially satisfied by fields, which depend 
only on either of the light cone coordinates z or z, for example the energy- 
momentum tensor and it's regularized powers ( in most of the models, we 
are going to look at, there are also other conserved currents ). If an infinite 
number of local integrals of motion survives the perturbation of the CFT, 
then we know that the massive theory is described by factorized S-matrices, 
of which some examples have been studied in chapter f|. There we made 
some assumptions, which resulted in the existence of an infinite number of 
conservedcharges P s , where s belonged to a certain subset of the natural 
numbers ^Hl. Suppose we perturb the critical action S* by a relevant scalar 

22 Recall that P- s is related to P+ s by parity, which we always assume to be a symmetry, 
allowing us to restrict ourselves to s > 0. 
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operator: 

S = S* -A J <f){z,z)d 2 z, (187) 

where the weight of <p is (h, h), so that the dimension of A is (1 — h, 1 — h). 
For this to be a relevant perturbation, we need y = 2(1 — h) > or h < 1. 
Charges, whose conservation is preserved by the perturbation, now become 

P s = j[T s+l dz +e,_itfc], (188) 

where T s and B s are local fields of spin s, satisfying the continuity equation 

dsT a+1 = d z Q a - 1 . (189) 

If S* described a free field theory, such a perturbation would be called 
super-renormalizable and adding the finite number of terms with dimensions 
smaller or equal than that of <fi one obtains a finite theory without UV- 
divergencies. In particular this doesn't change the structure of the UV fixed 
point and the fields continue to have the same dimensions at short distances, 
although the large distance behaviour of the perturbed theory is very much 
different. We will assume that this is also true for our case, where we perturb 
around a theory, which is not free. 

Thus S* describes an UV-finite ( renormalized ) theory, which contains 
<p as one of its operators and we assume, that we know all its correlation 
functions ( as well as those of all other operators ). Our primary aim is 
to check whether there exist currents J{z), whose conservation survives the 
perturbation, i.e. whether it satisfies equ.( |189| ). 

The perturbed correlation functions of a particular operator J(z) are now 
given by 

< J(z,z) . . . >=< J(z) . . . > s * +A J d 2 Zl < J(z)(j)(zx,z l ) ...> s * +0(X 2 ). 

(190) 

If this integral were finite, it would be independent of z. Therefore any 
^-dependence can come only from possible singular points z — > z\. In their 
neighborhood, we can use the short distance expansion ( SDE ) : 

J(z)(j){z h z j ) = r ~ : j^+A=Ar fo(*i>5j)> (191) 
k l~ z j\ 
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where A = 2h and Aj and A& are the scaling dimensions of J and <pk- Since 
only non-integrable singularities will contribute, this requires Aj + A — A& > 
2. In a unitary theory all dimensions are > and therefore only a finite 
number of operators (fik will contribute to first order in A. 

Let us look at the specific example of the energy-momentum tensor. In 
this case the SDE is 

T(z)(f)(z 1 ,zi) = , k v 2 0(zi,zi) + —^—Wizuzt). (192) 

\Z — Z\) Z — Z\ 

Since we only want to check whether the current T(z,z) is still conserved 
and to avoid bothering with infra-red singularities, let us only calculate the 
derivative of T(z,z). For this purpose we use the equation 

dz{z - 0~ m_1 = i-d?5^\z - 0- (193) 

ml 

This can easily be proved, either remembering that the logarithm is the 2- 
dimensional Green function Alnr = — 2ii5( 2 \r), i.e. 

d z (h ln[(z - z') (z - z')] = -2t™<5 (2) (z - z') 

and taking derivatives thereof, or proceeding as follows. Let us regulate 
the UV divergence by a cut-off a , inserting the step function H(x), which 
vanishes for x < 0. Then 

f ,2 ~ H[zz-a 2 ] 

J dzd * z m + l f (*>*) = 

r l2 zS(zz - a 2 ) ,, . r f(ae lv ,ae~ tlp ) ,^ niS 
Hence T(z,z) satisfies 

d- [ d 2 z x < T(z)ip(z,z) ...>= [ d 2 Zl dj < ( - — h —,(t>{z u z 1 ) + — *— 
j j Viz — ^1/ ^ — 

y rf 2 ^i (/i2tt9 2 (5 (2) (2 - Zl ) + 2tt5 (2) (^ - Zi)^) < <p( Zl , zi) . . . > . (195) 
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We therefore immediately get the conservation law for the energy- momentum 
tensor - as expected, since the energy-momentum tensor must remain con- 
served - : 

&zT + d z Q = 0, (196) 

where 

e = ir\(l-h)<f)(z,z). (197) 

We see, that the term in the SDE equ.( |192| ), relevant for the conservation of 
the current J(z, z) is picked out by : 

d I J(z,z) = \7r<f ^J( z )0( Zl ,^). (198) 

JCz 2,711 

Remembering the discussion on radial quantization leading to equ. (|159|) , we 
realize that the r.h.s. of this equation is a commutator and we get the 
following suggestive form : 

dsJ{z,z) = [J(z,z),H int (z)}, (199) 

where 



H, 



hit \ 



A J dz 1 (f)(zi, z) 



It remains to be seen, if the r.h.s. of equ. (|198|) can be expressed as d z of some 
operator. If yes, the conservation law will continue to hold to first order in 
A. 

To provide an example of how this can be checked, let us take the energy- 
momentum tensor and it's powers. Define A to be the irreducible Virasoro 
modul with highest weight h = 0, to which the energy-momentum tensor 
belongs. Introduce Virasoro generators as in equ.(161), but via an expansion 
around an arbitrary point ( : 

+ 00 T 

TW ^(c^ (200) 

Then, using the commutation relation 

[L_!,0(C,?)] = 3c0(C,£), (201) 

it immediately follows for the operator dz in equ.( |198|) - with J(z) replaced 
by T(z) - 

dzL^A = L_i<9jA. (202) 
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Now we streamline our algebra following reference . Introduce a set of 
operators D n , n = 0, ±1, ±2, ... as 

D n A(z,z) = l^^cm-zT^z), (203) 

i.e. D n projects out the term proportional to (( — z) n+1 in the SDE of (f) and 
A. The following equations show, why the D n are of good use to compute dz 

d- z = D . (204) 

Also 

D-n-j =/||(c-*r"~V(c^), n>0, 

and using the residue theorem for the (n+l)-th order pole : 

D- n -J = -.d n z (j)(z,z). (205) 

Since T(z) is generated from the unity / by applying operators L_ n , we need 
the commutation relations between L_ n and D m . First remember 

[L_n, 0(C, C)] = {(C - zT +1 d< + (n + l)h(C - zrmCC). (206) 
Now compute the commutator of [L n , D m ] : 

L n D m A(z,z) = f —L n (j)(C,z) = 

Jz Z711 

(C - z) m A(z) =fP- ([Ln, <f>((,z)] + <f>((,z)L n ) (C - z) m A(z). 
Jz Ln% 

Therfore 

[L n , D m ]A(z,z) =f^[L n ,<f>(C,mC-z) m A(z) = 

j ^ ((C - zr +m+1 d ( + (n + l)h(( - z) n+m ) <f>((,z)A(z) = 
{- ((1 -h)(n + l) + m)} jf ^0(C,^)(C - z) n+m A(z). 
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Hence we get the commutator 



[A, An] = -{(1 -h)(n + l) + m}D n+m . (207) 

A trivial application is 

dgT(z, z) = \D L_ 2 I = X(h- 1)D_ 3 J = X(h - l)L-i<j>(z, z) (208) 

reproducing equ.( |196|) . 

A less non-trivial calculation is to check the conservation of higher powers 
of the momentum. Let us define a regularized square of T T±(z) =: T 2 (z) : 
as 

T 4 (z) = (L_ 2 L_ 2 /)(z) = jf rfC(C - ^)^ 1 T(C)T(z). (209) 
Now check it's conservation : 

d-T± = \D L„ 2 L_ 2 I = \{h - 1){D_ 2 L_ 2 + L_ 2 D_ 2 )I = 

= X(h - l)(2A 2 A a + -^-L z _ x )<t>. (210) 

o 

Therefore, due to the first term above, the r.h.s. is in general not a derivative 
of z. This is to be expected, since the existence or not of this conservation 
law is a dynamical question. Let us take as perturbation one of the fields (f>i^ 
of the unitary models with c < 1 as an example. It is degenerate at level 3 
and repeating the steps that led to equ. (|172|) for the present case, we would 
get the following null- vector equation: 

I £-3 - + 7-r It. r-^ 3 1 ] 01 3 (z) = 0. (211) 

V h + 2) (h + l)(h + 2) > v > 

Hence the term containing L_ 2 in equ. ( |210| ) can be eliminated in favor of the 
derivative L_i : L_ 2 L_i = L_iL_ 2 — L_ 3 and using for L_ 3 the null-vector 
equation, we get 

dzT^z,z) = d z e 2 (z,z), (212) 

with 

h + 2\ 2 6{h + l) -iy^i,3 

Finally let us check that, similarly to the first order perturbation, the per- 
turbation expansion to n-th order, which usually contains an infinite number 
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of terms, is here also drastically truncated. The n-th order term will have 
the form 

„ n 

X n J d 2 Zl ...d 2 z n < J(z)l[ ( f>(z j ,z j )...>s* (213) 

We easily see , that the condition for non-integrable singularities is now 
Aj+nA-Afc > 2n or Aj — (2 — A)n - A fe > . Since A < 2, this condition 
will eventually be violated and the perturbation expansion has to stop with 
a finite number of terms. 



6.1 Counting Arguments 

In principle we may now take specific models and start looking for surviving 
conservation laws. It turns out 0, that at the expense of some formalism 
this job can be significantly simplified. Let us then present, what is called 
Zamoldochikov's counting argument. It will enable us to find out which of 
the P s are still conserved, at least for small s, without having to compute 
explicitly the term S in equ. (|189 | ) . Once we know the conservation laws, we 



may link up with chapter |2| and lift a candidate S'-matrix for the model in 
question. 

Let us talk about the energy-momentum tensor T, since T and its regular- 
ized powers will provide the conservation of the momentum and it's powers. 
The same type of reasoning also applies to other conserved currents. 

Recall that A was defined to be the irreducible Virasoro modul with 
highest weight h = 0, to which the energy-momentum tensor belongs. That 
is, A is the infinite-dimensional space spanned by all the fields of the form 
L_ ni L_ n2 . . . L_ nk /, where rij are positive integers or zero. A may be decom- 
posed as 

oo 

A = A s , (214) 

s=0 

where the fields belonging to A s satisfy Y%=i n i = s - From the commutation 
relations for the L n we easily see that the A s are eigenspaces of L : 

L A S = sA s . (215) 

Thus all fields belonging to A s have conformal weight (s, 0) and therefore 
dimension and spin equal to s. Besides this all these fields depend only on z 
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and are thus analytic, satisfying c^A = : they all give rise to integrals of 
motion. However for our counting argument it is important to exclude fields, 
which are derivatives of others. These would lead situations like d z d^T = 
dzR, which do not correspond to conservation laws at all. They are contained 
in L_iA. Therefore let us define a new space A, where these fields are divided 
out : A = A/L_]A. This space also has the decomposition 

oo 

A = A s , L A S = sk s . (216) 

s=0 

We can now take advantage of character formulas to obtain the dimensions 
of A s : 

oo 

q s dvmk s = (1 - q) Xo (q) + q, (217) 

s=0 

where Xo(q) is the character of A defined as 

oo 

Xo(q) = E q S dimk s . (218) 
Equ.(^T^) can easily be shown as follows. First note that 
dim ((L_iA) s ) = 



dimk s -i , for s > 1 
s = 1, 



since L_\I = 0. Therefore : 



q s dim ((L_iA) s ) = V] q s dimA s _ 1 = q^^q s 1 dimk s _i = 

s=0 s=2 s=2 

oo 

?E? S dimA s - q = qxo(q) ~ Q- 

s=0 

and equ.( |217D follows. If c doesn't belong to the degenerate set given by 
equ.( |175|) or equ.( |177| ), the dimension of A s , according to its definition, equals 
p(s) - the number of partitions of the integer s. Equ. (|182p gives then 

oo oo 

J2q s p(s) = l[(l-q n r 1 = Xo(q), (219) 
s=0 n=l 
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so that finally we get t^l: 

oo 

q s dim(A s ) = 1 + q 2 + q 4 + 2q 6 + . . . . (220) 

s=0 

Now we have to ask, which of these conserved currents survive perturba- 
tion ? Suppose we add a relevant perturbation \<f>(z,z). (f>(z,z) is the highest 
weight vector of the irreducible modul $0$, generated by all the fields of 
the form L_ ni . . . L_ nfe L_ mi . . . L_ m ,0, where rij, rrij are positive integers or 
zero. As A, $ ( and $ ) may be decomposed as 

oo 

$ = $ s , L $ s = (A + s)$„ X $ s = A$ s . (221) 

Now the fields e A s do not satisfy d^T^ = 0, but 

q_ t ( k) = x^} 1 + 0(X). (222) 

Here are local fields belonging to ^spi and for simplicity assume, that 
only the first order term in A contributes 1^1. We use this equation to define 
the linear operator 

% : A s -> $ s _i. (223) 

Let us consider the space $, where we factored out the derivatives : $ s = 
$ S /L_i$ s _i. Consider now the mapping .M s from A s — > $ s _i. A4 S is im- 
plemented by Al s = H S D ^ S , where Il s is the projector from $ s to $ s and 
Dq }S is <9j restricted to A s . Since in $ we factored out the derivatives, all the 
fields satisfying 

d^T s+1 = d x e a - h (224) 

are mapped into the null element of i.e. all elements T s+1 satisfying 

equ.( [224| ) belong to the kernel of A4 S . If dim A s+1 > dim $ s , then Ker M. s ^ 
and we have conserved charges surviving the perturbation. 

The dimension of $ s can be computed as that of A s and we obtain: 

oo 

]T q A+s dim$ s = (1 - q)xA(q), (225) 

s=0 



23 For the degenerate cases, we have to subtract invariant subspaces and use equ.( |186| ). 
Even then equ.(217) is still valid for s < m(m — 1) 



21 



See ref.Q for a discussion when this is not true. 
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where Xa{q) is the character of the modul with highest weight (h, h) 
Let us apply this scheme to specific cases. 

1) Lee- Yang edge singularity and the associated non- unitary 
series 

The Lee- Yang singularity describes the critical behavior of an Ising model 
in a purely imaginary magnetic field ih. Above a critical value h c the zeroes 
of the partition function condense on the imaginary h-axis and the following 
Ginzburg-Landau lagrangian describes this behavior: 

C = j d 2 x[^{d^{x)f - i(h - h c )<p{x) - i<p{x) 3 ). (226) 

This is clearly a non- unitary theory and at the critical point h = h c <p(x) 
is the only relevant operator. In ref.JJ it was shown that there exists only 
one model with these properties and no extra symmetries. With the aid of 
equs. (|175|) and equ. (|176|) one easily identifies it to have p = 2 and p' = 5 
with central charge is c = c(2, 5) = —22/5. It contains only two primary 
fields : the identity and the field ip = ^(1,2) = ^(1,3) with weights (0, 0) and 
(—1/5,-1/5) respectively. The negative values reflect lack of unitarity here. 
We now check for conservation laws computing the dimensions of the 



spaces A s and $ s , using equ. fl217j ) and equ. (|225|) . One finds, that the dimen- 



sion of A s exceeds the one of $ s by one unit for 

s = l,5,7. (227) 

These are the first 3 conservation laws found in section |2|. This reasoning has 
been extended^] to the whole series of non-unitary models p — 2,p' — 2N + 3 
with 

2M6Y + 5) , x 

c(2,2iV + 3 = ^ 228 

v ' ; 2Y + 3 K J 

perturbed by the field ^(1,3) • 
2) Ising model. 

It's conformal anomaly is c = c{m = 3) = 1/2. This model contains 

three spinless primary fields : / = <^(i,i),cr = ^(1,2) and e = ^(1,3) with 



2r, 



Notice that the identity is the only primary field satisfying L_i/ = 0, which accounts 



for the difference between equs.(217) and (225) 
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conformal weights (0, 0), (1/16, 1/16) and (1/2, 1/2), which are identified with 
the identity operator, the spin density (magnetization) and the energy density 
respectively. All the other local fields in this model are obtained by applyig 
the left and right Virasoro generators L n and L n with n < to these primary 
fields. 

Suppose we perturb this model by a magnetic field, which breaks the 
Z 2 invariance and couples to the magnetization a(x) = (pn^x). The total 
action is 

S[)f = S* + h J d 2 xa(x). (229) 

The dimensions of the spaces A s and $ s can be calculated, using equ.( pl7|) 
and equ.( [225p . For the characters, we use equ.( |186| ). The relevant formulas 
are 



1 [ OO OO 1 

Xo(q) = XiM = o II (! + <1 H+1/2 ) + II (1 - ( l n+1 ' 2 ) > (230) 

Z ln=l n=l J 



OO OO 

Xi/m{q) = XiM = IK 1 " f n+1 T l = 1 1/W II (1 + 1 U )- (231) 

n=0 n=l 

Computing now the dimensions, we check that the dimension of A s exceeds 
the one of $ s by one unit for 

5 = 1,7,11,13,19. (232) 

For larger values of s the dimension of $ s is greater or equal and nothing 
can be concluded. However we obtained five nontrivial conservation laws 
for the Ising model in a magnetic field and one conjectures the infinite set, 
corresponding to all integers s, relatively prime to 30 as found in section 0. 



3) Z(N)-models. 

Using again the by now familiar scheme let usconsider the Z(3) model 
perturbed by the field Wi,2); which is the most relevant thermal operator. 
We get conservation laws for s = 1, 2, 4, 5, 7, 8, which are the first six integers 
of the infinite series s ^ 0,mod(3). 

This is a particular case of the following general results obtainable using 
the machinery developped in the lectures. 
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The UV limiting CFT of the Z(N) models exhibit primary fields of integer 
spin s — 3, 4, . . . , N, which depend only on the variable z - and analogously 
for ~z. These are therefore conserved currents and the Virasoso algebra is 
enlarged to a W(N) algebra, containing the energy-momentum tensor and 
those currents 0. This additional symmetry is also suggested by the coset 
construction |J, from which the Z(N) models may be obtained, using the 
A N _i Lie algebra. The central charge of the W(iV)-models depends on two 
integers: 

e(p,A0 = (A-l)(l-^±i>). (233) 

For p = N + 1 we obtain our Z(N) series. The generalization of the 
f(i,2) £ Z(3) is a primary field of weight ( h 2 +2 ? ) ' w ^ich preserves the 
global symmetry and is the most relevant thermal operator. Here h A = N 
is the so-called dual Coxeter number of the Lie group -Ajy-i- Perturbing 
the UV limit with this operator, we obtain our massive Z(N) models with 
conservation laws P s , where s takes the exponents of Ajv-i ( which are all 
integers 1,2, ... ,N ), repeated modulo the Coxeter number h? , i.e. modulo 
N. 
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